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ABSTRACT 

The gamma-ray binary LS I +61 303 is composed of a Be star and a compact companion orbiting in an eccentric orbit. Variable flux 
modulated with the orbital period of ~26.5 d has been detected from radio to very high-energy gamma rays. In addition, the system 
presents a superorbital variability of the phase and amplitude of the radio outbursts with a period of ~4.6 yr. We present optical 
photometric observations of LS I +61 303 spanning ~1.5 yr and contemporaneous Ho: equivalent width (EWyia) data. The optical 
photometry shows, for the first time, that the known orbital modulation suffers a positive orbital phase shift and an increase in flux for 
data obtained 1-yr apart. This behavior is similar to that already known at radio wavelengths, indicating that the optical flux follows 
the superorbital variability as well. The orbital modulation of the EWy{ a presents the already known superorbital flux variability but 
shows, also for the first time, a positive orbital phase shift. In addition, the optical photometry exhibits a lag of ~0.1-0.2 in orbital 
phase with respect to the EWu a measurements at similar superorbital phases, and presents a lag of ~0.1 and ~0.3 orbital phases with 
respect noncontemperaneous radio and X-ray outbursts, respectively. The phase shifts detected in the orbital modulation of thermal 
indicators, such as the optical flux and the EW\{ a , are in line with the observed behavior for nonthermal indicators, such as X-ray or 
radio emission. This shows that there is a strong coupling between the thermal and nonthermal emission processes in the gamma-ray 
binary LS I +61 303. The orbital phase lag between the optical flux and the EWu a is naturally explained considering different emitting 
regions in the circumstellar disk, whereas the secular evolution might be caused by the presence of a moving one-armed spiral density 
wave in the disk. 
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1. Introduction 


LS I +61 30 3 is one of the five gamma-ray binaries currently 
known fe.g.. lDubusll2013| : IParedes et alJl2013h , It is composed 
of an optical star (V~10.7) with spectral type B0 Ve, therefore 
presenting a circumstellar disk, and a compact companion that is 
probably a p ulsar orbiting in an eccentric orbit with e in the range 
0.54-72 (seelHutchings & Cramntonll 1 981b IParedes & Fig ueras 


1986 : ICasares et alj 120051 : iDhawan et alj 120061: Aragona et al 
2009). The distance to the s ource is estimated to b e ~2 kpc 
based on Hi measurements (Frail & Hiellming 1991). The bi¬ 
nary _system has an orbital period of 26.4960 ± 0.0028 d 
dGregorvl2002lL and variable flux modulated with the orbital pe¬ 
riod has b een detected from radio to very high-energy gamma 
rays (e.g., Taylor&Grggory ^982; Mendelso n & MazebI 1989: 


Paredes et alJ I 19941 1 19971: lAbdo et alJl2009l: lAlbert et aljl 2009). 


In addition to the orbital modulation, LS I +61 303 also exhibits 
a superor bital modulati on of radio outbursts in ~4,6 yr (first de- 
tected bv lParedesI 19871 and lGregorv et al.l 1989b . iGregorvI (12002b 
found a superorbital period of 1667 + 8 d in the amplitude and 


* Tables 1 and 2 are only available in electronic form at the CDS 
via anonymous ftp to cdsarc.u-strasbg. fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/ 

** Serra Hunter Fellow. 


in the orbital phase of the radio outbursts, leading to a drift of 
the outburst maxima from orbital phase ~0.4 to ~0.9. For ref- 
erenc e, the periastron passage takes place at o rbital phase 0.23- 
0.28 (ICasares e t al. 2005: lAragona et akll2009b . The zero orbital 
and superorbital phases are set at ID 2,443,366.775 ( Gregor y 

!^ol . 


IZamanov et ali(l2013l) have recently detected the ~4.6 yr su¬ 
perorbital period in the equivalent width of the H a e miss ion 
line, confirming the earlier evidences reported in IZamanov et al.l 
(11999b . This suggests that the superorbital variability is related to 
periodic changes in the mass-loss rate of the Be star and/or varia¬ 
tions in the circumstellar disk (see lMassi & Torricelli-Ciamponil 
20141 a nd references therein for alternative interpretations). 


Li et al.l (12012 ) discovered the superorbital modulation of 


LS I +61 303 in X-ra ys (3- 30 keV), with a period compatible 
with that found by Gregory d2002b from radio measurements. 


__ . 

IChernvakova et all (12012b also reported the superorbital modu¬ 
lation of the X-ray emission, and using contemporaneous X-ray 
and radio observations showed that radio outbursts lag the X-ray 
outbu rsts by ~0.2 in orbital p hase along the entire superorbital 
cycle. [Ackermann et al.l (12013b found a sinusoidal variability in 
the >100 MeV y-ray flux consistent wi th the radio superorbital 
period. Finally, [Zaitseva & Borisov ! (2 003ft found long-term vari- 
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Fig. 1. Optical photometric lightcurve (top, red) and EWn a (bottom, blue) of LS I +61 303 plotted as a function of MJD (left) and folded with 
the orbital phase (right). Error bars represent l-cr uncertainties. The solid curves in the right-hand panels represent sinusoidal fits to the orbital 
variability, whose parameters can be found in Table [3] Two cycles are displayed for clarity. 


ability in the mean V -hand magnitude of 0.07 mag, but no search 
for superorbital variability was conducted. 

In this work, we present optical photometric observations of 
LS I +61 303 spanning ~1.5 yr, and contemporaneous Her equiv¬ 
alent width (EWfta) observations. We report the discovery of an 
orbital phase shift and a variation in the orbitally modulated flux, 
very similar to the well-known superorbital behavior of the or¬ 
bital modulation in radio. In addition, we also detect a phase 
shift in the orbital modulation of the EWu a , to our knowledge 
reported here for the first time. We discuss the obtained results 
in the context of the radio and X-ray superorbital variability, and 
conclude that there is a strong coupling between the thermal and 
nonthermal emission in LS I +61 303. 


2. Observations and data reduction 

We performed optical photometric observations of LS 1+61 303 
with the robotic Telescope Fabra-ROA Montsec (TFRM; see 
iFors et akll2013l) . The TFRM is installed at the Observatori As¬ 
tronomic del Montsec (Lleida, Spain). The main specifications 
are: corrector plate of 0.5 m aperture and 0.78 m primary mir¬ 
ror, refurbished Baker-Nunn Camera for routine CCD robotic 
observations, focal ratio f/0.96, 4.4° x 4.4° field of view with a 
pixel scale of 3.9"/pixel, passband filter SCHOTT GG475 (A > 
475 nm.), and custom CCD based on FLI ProLine 16803 with 
quantum efficiency of 60% at 550 nm (IFors et al.ll2013l) . 

The observations span from 2012 July 31 to 2014 March 7 
(2 seasons) with 71 nights of good datajpreliminary results on 
season 1 data were presented in lParedes-Fortunv et al.l20l4b . We 
observed the target around 20 times per night with exposures of 
5-10 s, as a compromise between good signal to noise ratio and 
avoiding the nonlinear CCD regime. 

We conducted the data reduction and analysis using a 
pipeline developed in Python following these steps: standard 
calibration of the images using IRAlQ (implemented through 
PyRAF0), aperture photometry using the PHOT package from 
IRAF (PyRAF) with an aperture radius of 7.5 pixels, and cor- 

1 IRAF is distributed by NOAO, which is operated by AURA, under 
cooperative agreement with NSF. 

2 PyRAF is a product of the Space Telescope Science Institute, which 
is operated by AURA for NASA. 


rection of the lightcurves using a weighted _average differen¬ 
tial magnitude correction method based on lBroeg et al.l (l2005i) . 
The lightcurve of LS I +61 303 has been corrected using 145 
weighted reference stars up to 4 mag fainter than the target and 
as close as possible to it (< 0.25°). Finally, we averaged the cor¬ 
rected magnitudes of the target on a nightly basis, and we used 
an artificial offset to a mean magnitude of 10.7 mag. The nightly 
photometric uncertainties are estimated as the standard devia¬ 
tion of the magnitudes of the target obtained from the individual 
images for each night. The obtained values are typically in the 
5-10 mmag range. 

We also considered contemporaneous EW^a data of 
LS I +61 303 obtained using FRODOspec on the robotic 2.0- 
m Liverpool telescope at the Observatorio del Roque de Los 
Muchachos (La Palma, Spain). The data span from 2012 July 
2 to 2014 January 15 with 104 measurements, with one 600 s 
exposure per night. The data have been obtained with the same 
setup, and reduced in the same way, as in lCasares et al.l d2012l) . 
We assumed that the uncertainties of the EWu a are at the 10% 
level. The adopted convention is positive EW for emission. 

The data (Tables 1 and 2) are available online at the CDS. 

3. Results 

The optical photometric lightcurve and EW\r a of LS 1+61 303 as 
a function of MJD and folded with the orbital phase are shown in 
Fig. ID As can be seen, the data cover two different seasons (cen¬ 
tered in autumn) and show the already known orbital variabil¬ 
ity (see, e.g., IZamanov et al.l (20141) . The folded lightcurves for 
the first and second season and their sinusoidal fits are shown in 
Fig. [2] The fitted sinusoidal parameters (amplitude, zero phase, 
and mean value) for the different data sets, considering the data 
uncertainties, are quoted in Table [3] The l-cr uncertanties in the 
fitted parameters quoted in Table [3] are estimated from the co- 
variance matrix. 

The fits to the orbital variability of the optical photometry re¬ 
veal, for the first time, a positive orbital phase shift of 0.16+0.03 
at 5.2 cr confidence level (c.l.) and an increase in the average op¬ 
tical flux of 15 + 4 mmag at 4.2 cr c.l. between the two obser¬ 
vational seasons 1-yr apart. These trends are reminiscent of the 
superorbital trends found at radio and X-ray wavelengths, and 
indicate that the optical flux follows the superorbital variability 
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Fig. 2. Optical photometric lightcurve (left, red) and EW HtK (right, blue) of LS 1+61 303 folded with the orbital phase for the first (top) and second 
(center) seasons, and their sinusoidal fits to the orbital variability (bottom; see Table 0. Error bars represent 1-cr uncertainties. Two cycles are 
displayed for clarity. The shifts in the orbital modulations between both seasons are clear. 


Table 3. Fitted sinusoidal parameters to the orbital variability of the optical photometry and the EWn a of LS I +61 303 for the total data set 
(Fig. 03 and for the two different observational seasons (Fig. Q. The parameters are those of the functions: optical m = -A cos (In [ip - (J>q\) + C 
and EW\\ a = A cos (2 n \tj> - 0o]) + C. Phases were computed using an orbital period of 26.4960 d and phase zero at JD 2,443,366.775. 


Data set 

MJD optical 

A optlcal (mmag) 

.optical 

C°P t,cal (mag) 

MSD l:w "“ 

a ew H" (A) 

<P, 

c EW »" (A) 

Total 

56139-56723 

28 + 3 

0.73 ± 0.02 

10.701 + 0.002 

56110-56672 

1.2 ±0.2 

0.56 ± 0.03 

11.5+0.2 

Season 1 

56139-56335 

29 + 3 

0.68 ± 0.02 

10.704 ± 0.003 

56110-56298 

1.4 ±0.3 

0.54 ± 0.03 

11.2 + 0.2 

Season 2 

56594-56723 

24 + 4 

0.85 + 0.02 

10.689 ± 0.003 

56603-56672 

1.3 ±0.3 

0.76 ± 0.05 

13.0 + 0.3 


as well. We note that the folded lightcurves of the two individ¬ 
ual seasons (see Fig. (2]-left) show less scatter than the folded 
lightcurve including all the data (see Fig.QJ. which span ~0.35 
superorbital phases. 

Similarly, the fits to the orbital variability of EW\\ a present 
for the first time a positive orbital phase shift of 0.22 ± 0.05 
at 4.2 cr c.l. Further, we confirm the previously claimed su¬ 
perorbital modulation of the mean EW\\ n , with a variation of 
1.8 + 0.3 A at 5.7 cr c.l., between our two observational sea¬ 
son^. In addition, the orbital modulation of the optical photom¬ 
etry exhibits a lag of ~0.1-0.2 in orbital phase with respect to the 
orbital modulation of E VLh,, at similar superorbital phases (see 
Fig. [2] and Table[3}. 

To better display the superorbital variability, color maps of 
the optical photometry and EW\\ a as a function of the orbital 
phase and superorbital cycle are shown in Fig. [3] The dotted 
red and blue lines represent the orbital phase drifts of the cor¬ 
responding maxima along the superorbital cycle for the contem¬ 
poraneous optical photometry and EWn a presented here (super¬ 
orbital cycle 7). The contemporaneous radio and X-ray fluxes 
of the previous superorbital cycle (6) from IChernvakova et all 
( 2012 ) are shown as black and green dotted lines. The optical 
photometric observations show a lag of ~0.1 in orbital phase 
with respect to the radio outburst, and ~0.3 with respect to the X- 

3 Despite the re lative ly poor sampling at equivalent superorbital 
phases, the lZamanov et al] (2(113 ) data reveal a compatible behavior. 


ray outbursts for equivalent superorbital phases one cycle apart. 
The EWhiy maxima occur at orbital phases similar to the radio 
outbursts for similar su perorbital phases one cycle apart, while 
in the averaged data of lZamanov et alJ (2013. 2014 ) the maxima 
occur during the rising of the radio flux density. The different 
slopes in Fig.[3]might suggest a change in the behavior of the su¬ 
perorbital modulation between different superorbital cycles, al¬ 
though we caution that the results are still compatible at 2 cr c.l. 


4. Discussion 

The observations of LS I +61 303 reported here show, for the 
first time, orbital phase shifts in the maxima of the orbitally 
modulated optical flux and EW^a for different phases of the su¬ 
perorbital cycle (see Figs. [2] and [3] and Table [5]). These results 
extend previously reported links between the Her emission and 
multiwavelength properties but obtained from the average or¬ 
bital modulation, i.e., without considering the superorbital vari¬ 
ability (see, e.g.. iZamanov et akll2014l) . The results reported here 
are now discussed in the context of the known multiwavelength 
behavior of the orbital and superorbital variability of the source. 

The orbital variability of the gamma-ray binary LS I +61 303 
is wavelength dependent. After periastron at phase 0.23-0.28, 
there is first a maximum in X-rays (good tracer of the nonther¬ 
mal emission), followed by a maximum in EWn a (tracer of the 
outer disk conditions) and in radio emission (produced mainly 
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Fig. 3. Color maps of the optical photometry (left) and EWa a (right) from LS I +61 303 as a function of the orbital phase and superorbital cycle. 
The red crosses and blue stars correspond to the phases of the maxima of the sinusoidal fits to the orbital variability of the optical photometry and 
EWu a , respectively, for the first and second seasons (see Table0. The dotted color lines represent the orbital phase drift of the emission peaks 
caused by the superorbitaLvariability for the optical photometry (red), EWu a (blue), radio (black), and X-rays (green). The radio and X-ray drifts 
are taken from Fig. 3 of tChernvakova et~al . ( 2012 1 using data of the previous superorbital cycle. 


outside the binary system). Finally, a maximum in optical flux 
is observed, with a 65% contribution from the Be star plus 35% 
from the inner circumstellar disk (iCasares et al .1120051) . The cir- 
cumstellar disk is likely to be perturbed and/or partially dis¬ 
rupted by tidal forces and the putative pulsar wind ram pres¬ 
sure. These would trigger significant changes in the structure of 
the Be disk, especially around periastron passage, reducing its 
emitting area and total optical emission. This is approximately 
what is observed in the data from season 1, but the orbital phase 
shift in season 2 clearly shows that there is more than purely 
orbit-induced variability. In fact, the superorbital variability ob¬ 
served in EW\r a has been associated with periodic changes on the 
Be star envelope and its circumstellar disk (e.g.. tZamanov et aki 
119991) . which would trigger the superorbital variability observed 
at other wavelenghts. If this is the case, since the orbital mod¬ 
ulation seen in X-rays or radio suffers a phase drift along the 
superorbital cycle, we should detect the same effect in optical 
flux and EWu a , as reported here for the first time. Thus, it is 
clear that there is a strong empirical coupling between the ther¬ 
mal (optical) and the nonthermal (X-ray and radio) emission in 
the gamma-ray binary LS I +61 303 at both orbital and superor¬ 
bital scales. 

The orbitally modulated optical flux and EW\\ a present dif¬ 
ferent behavior between them. The optical flux shows ~0.06 mag 
modulation, representing ~6% in total flux or ~16% in disk flux 
(and probably projected area), while EW\\ a shows ~30% vari¬ 
ability. This implies that external parts of the disk are more per¬ 
turbed or disrupted than the inner parts of the disk, as one would 
expect if these perturbations (or disruptions) are due to the influ¬ 
ence of the compact object as it approaches periastron (where 
the perturbations are caused by tidal forces and/or ram pres¬ 
sure). In addition, the ~0.1-0.2 phase lag implies that the ex¬ 
ternal parts are perturbed or disrupted before the internal parts, 
but they recover earlier as well, probably building up from ma¬ 
terial shocked before periastron passage. The inner disk would 
only recover close to, or even after, apastron. The secular evo¬ 
lution of this behavior along the superorbital cycle could be due 
to the presence of a moving one-armed spiral density wave in 


the disk, as suggested bv lNegueruela et al. ( 1998 ) for traditional 
Be/X-ray binaries. 
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